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Motivated by recent experiments on DNA torsion-force-extension characteristics we consider the 
writhing geometry of open stiff molecules. We exhibit a cyclic motion which allows arbitrarily large 
twisting of the end of a molecule via an activated process. This process is suppressed for forces larger 
than femto-Newtons which allows us to show that experiments are sensitive to a generalization of the 
Calugaxeanu- White formula for the writhe. Using numerical methods we compare this formulation 
of the writhe with recent analytic calculations. 



I. INTRODUCTION 

Recent experiments in which DNA molecules are ma- 
nipulated 0, 0, El with the help of magnetic beads 
have led to a renewed interest in the statistical me- 
chanical properties of the torsionally stiff wormlike chain 
lEIEIEEEIS- I n the experiments a bead is attached to 
the end of a long DNA molecule; the other end remains 
stuck to a surface. The bead is held in a magnetic trap 
which allows the simultaneous application of a force, F 
and couple, T; one measures the mean distance of the 
bead from the surface as a function of F and T. In an 
ingenious theoretical paper it was shown 0] that the sta- 
tistical mechanics of this problem are related to the prob- 
lem of a quantum particle in a magnetic field. However, 
a crucial assumption was made in the formulation of the 
writhing geometry of the polymer. 

The concept of writhe was originally defined in the 
context of studies on closed ribbons where it forms one 
part of a topological invariant, the linking number [sLllOj. 
In this paper we show how the writhe can be usefully 
generalized to study the geometry (and the mechanical 
response functions) of open DNA molecules. This gener- 
alization introduces end corrections to the Calugareanu- 
White formula for the writhe. We then perform simula- 
tions for the writhe distribution of an open chain which 
we compare with the analytic theories. We show that 
experiments involving manipulation of DNA with beads 
in unconfined geometries have unbounded fluctuations in 
the measured torsional angle. Under an exterior torque 
a bead can rotate an arbitrarily large angle. When the 
DNA is under tension this rotation is an activated pro- 
cess with jumps of 47r in the mean angle. 

Contrary to the calculation performed by Bouchiat 
and Mezard |s| we find that this unbounded response is 
not removed by the introduction of an intermediate cut 
off. However applications of tensions larger than femto- 
Newtons suffice to render the problem finite in practice 
PH. We find that their regularized expressions strongly 
overestimate the magnitude of writhe fluctuations for 
tensions larger than femto-Newtons. 

In order to interpret experiments in which the molecule 
is under strong tension Moroz and Nelson used a Mongc 
representation to perform a calculation of writhe fluctua- 
tions . We numerically investigate the validity of their 



results. We show that a strong correction to scaling can 
be expected due to the formation of rare loops, which give 
however exceptionally large contributions to the writhe. 
We find that the window of forces in which simple an- 
alytic theories, based on the Monge representation, are 
valid is rather small. 

In single molecule experiments, self-avoidance confines 
the polymer to a single invariant knot. Analytic theories 
are unable to estimate the error in writhe due to sum- 
ming over both knotted and unknotted configurations. 
We investigate this question numerically. 



II. WRITHE AND LINKING NUMBER 



Closed curves 



The linking number Ck of a closed ribbon is an integer 
topological invariant [E El • It can be decomposed into 
two parts the twist, Tw and the writhe, Wr: 



Ck = Wr + Tw. 



(1) 



This decomposition is useful because the writhe is a func- 
tion of the centerline r(s) of the ribbon. 



Wr = i- / </« 

An 



ds' 



r(s) — r(s')| : 



drjs) x dr(s') 
ds ds' 



(2) 



The linking number is invariant under deformations of 
the shape which do not introduce self intersections. If 
during a deformation the centerline r(s) crosses itself 
there is a discontinuity of 2 in Ck and thus, as Tw is 
continuous, a discontinuity of 2 in Wr. Fuller 0, Iu3 ] 
showed the integral of equation (J2J) could be simplified 
and introduced the expression 



Wr F = ±- 
2tt 



e 2 • (t x t) 
1 +t • e. 



ds 



(3) 



where e z is the direction at both extremities of the open 
chain. In this simplification information is lost so that 
Wr and Wr F are related by the equation 



Wr - Wr b = mod 2 



(4) 
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For notational convenience let us introduced the angles, 
Xc = — 27rWr and xf — —2TrWr F . A much more di- 
rect approach to Fuller's result is possible 0] by noting 
that the writhe of a stiff polymer is closely related to the 
geometric anholonomies discovered by Berry |15| in wave 
phenomena. 

Equation J2J) has a simple geometric interpretation^^]. 
The projected path of a chain on a plane with normal u 
can intersect itself. Each crossing is assigned a number, 
±1, according to the handedness of the intersection. The 
sum of these numbers n(u) is the writhe number for direc- 
tion u and the expression, eq. J3J), is equal to the average 
of n(u) over all directions. 

B. Bead rotation is given by an extended definition 
of writhe 

The above definition of the writhe, with its emphasis 
on it being part of a topological invariant hides, to some 
degree, the interpretation of xc as an angle of rotation in 
many experimental situations. We shall now show that 
an extended definition of the writhe, based on equation 
(J2J) is the writhe contribution to the bead rotation that 
is measured in DNA twisting experiments. Rather simi- 
lar arguments have also been given in Q for a polymer 
between two planes, we give hear an extended derivation 
to point out the exact limits of the result. 

Consider the planar ribbon in figure ^ The linking 
number, Ck, is zero thus 

Wr + Tw = 0. (5) 

We shall use this closed ribbon as a reference configu- 
ration to calculate the writhe of an open filament for 
which the initial and final tangents are parallel i.e. for 
t(0) = t(L), where L is the polymer length. We do this 
using the geometry of figure The polymer V which is 
also a ribbon is embedded in a construction consisting of 
two long straight sections Si, S2 and a closing loop C. We 
shall apply eq. @ to the constructions of figures ^and [21 
then take the length of the straight sections to infinity. 

There are beads attached to each end of the polymer, 
joining on to the straight sections. These two beads are 
our experimental reference. We shall hold the lower bead 
Bi stationary and let the upper bead B2 rotate due to the 
writhe of the polymer. The final part of our construction 
is the "twist absorber" T. This is imagined as being a 
joint, or section of the straight ribbon which twists freely 
to absorb the writhe generated by the polymer. Let the 
twist of the polymer be zero so that all twist appears in 
T. 

Start from the reference state of figure ^ an d deform 
continuously to an arbitrary state figure|21without gener- 
ating a self intersection of the construct. The writhe can 
be calculated with the classic double-integral. Clearly 
since Wr + Tw = the writhe which is generated goes 
into twisting the region T of the chain. The total twist- 
ing angle is just — 2nWr = xc- Now let the straight 




FIG. 1: A flat ribbon with zero linking number. 

sections go to infinity. In this limit the contribution to 
the integral of the section C vanishes. Thus the experi- 
mental rotation angle is given by the xc for the extended 
construct with the two straight line sections extending to 
infinity |3lj| . 

Consider generating an ensemble of chains with some 
arbitrary algorithm. The distribution of the configura- 
tions of the chains is independent of the dynamic process 
creating them if they are subject to Boltzmann statistics. 
If we now choose unknotted configurations and demand 
that they be created via a process which preserves the 
linking number we conclude that we can uniquely deter- 
mine the rotation angle from the writhe of the extended 
construct. A crucial part of this argument is the restric- 
tion of the construct of figure to the sector Ck = 0. 
We discuss the importance of this restriction in the next 
section. 



C. Choice of linking sector 

In the deformation from figure ^ to figure |3 we wish 
to conserve the linking number Ck. To do so we must 
firstly forbid self crossing of the polymer configuration. 
In any experiment with DNA in the absence of topology 
changing enzymes this restriction in reasonable. However 
this condition is insufficient to conserve the topology of 
the entire construction. We must also forbid the cross- 
ing of the real polymer with the imaginary line from the 
end of the chain to infinity. In a recent numerical paper 
Q this was achieved by grafting the ends of the polymer 
to an external surface rather than beads. Experimen- 
tally we would stop this crossing by introducing a steric 
hindrance: A fine fiber in the neighborhood of the manip- 
ulating bead would work perfectly. We show below, how- 
ever, that even this level of steric hindrance is not needed 
in the experiments as they are presently performed. 

If the chain is allowed to bend back in such a way that 
is passes through the line to infinity there is a discon- 
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FIG. 2: Deform the ribbon of figure^as in the diagram: The 
initial and final tangents of the polymer section V are parallel. 
Long arms S\ and 52 attach the polymer to a closing loop 
C. The experimentally important reference beads B\ and B2 
are attached to the two ends of the polymer. The twist is 
absorbed exclusively in the section T . Note the beads are not 
in general aligned vertically. 




FIG. 3: Demonstration of a cyclic motion of DNA leading to 
introduction of a twist of 47r. In the first and fourth figure the 
shape of the molecule is identical. Between the second and 
third figure the vertical extension of the bead passes through 
the molecule leading to a discontinuity in the linking number. 
The tricky three dimensional geometry is most easily under- 
stood by repeating this sequence with the help of a belt or 
ribbon. 



tinuity of ACk = ±2. This allows one to construct a 
cycle of motions in which the shape of the filament un- 
dergoes a cycle coming exactly back to its original shape 
while the bead undergoes a rotation of Att about the ver- 
tical axis, figure [3J This result is the origin of a common 
amphitheater demonstration of the importance of spinor 
representations of the rotation group: If one holds a plate 
horizontally in the palm of one's hand one can spin it 
about a vertical axis by performing a suitable contortion 
of the arm. Against all intuition the plate can be turned 
an arbitrarily large angle. Photographs of Feynman per- 
forming this demonstration are to be found in |17| . Each 
cycle of the arm again gives rise to a rotation of An of 
the plate. 

This unbounded torsional fluctuation is an activated 
process: When the DNA is under tension the bead moves 



against the magnetic force a distance comparable to its 
diameter in order to force a loop over the point of at- 
tachment. For micron sized beads this introduces a char- 
acteristic scale of the force of fcsT//im ~ fN. For forces 
of more than a few femto-Newtons the process is expo- 
nentially rare and suppressed. The natural time scale for 
attempts at crossing the barrier is given by the Zimm 
time of the bead t ~ cPrj/ksT with d the bead size, 
which is seconds. 

This force is extremely low, the widely studied 
crossover from random coil to semiflexible behavior in 
the force-length characteristics of DNA occur at a force 
scale of ksT/ip with i v « binm the persistence length 
of DNA. This intrinsic force scale for DNA much higher 
than that which is needed to conserve linking number 
with a 1/Ltm bead. We conclude that most experiments 
with tense DNA are performed in a regime where a high 
energy barrier leads to conservation of the linking num- 
ber of the construct of figure 2. However, there is always 
a small probability of passage by this barrier so that un- 
der torsion the true steady state of a DNA molecule is 
a state of cyclic motion in which the bead rotates by a 
series of activated jumps. 

Because of these considerations we assume that all 
topological invariants are conserved during experiments, 
including linking number and the knot configuration of 
the extended construct of figure 2. 

D. Analytic expressions for the extended writhe 

From the above discussion we are lead to the calcu- 
lation of the writhe of the extended construct, figure 2. 
The double integral, eq. J2J, on the interior of the chain 
is identical to that commonly used for closed DNA, we 
denote the corresponding angle by Xint ■ 

The integral from one point on the linear extension to 
infinity plus the second point on the polymer can be sim- 
plified as follows: The integral is evaluated by noticing 
that it is the spherical area, swept out by the vector 



as s and s' vary. Consider, now, a chain, r(s). Place 
the origin at r(0). Let us now place s in the interior of 
the polymer and s' on the extension of the polymer to 
infinity which is directed in the direction e z . Consider, 
now, s' = 0. The curve u(s) = e(s,0) is a spherical 
curve. As we now let s' vary from to infinity we sweep 
out the area between u and the point e 2 . This spherical 
area can be written as 

f u x du/ds 
Xend = e z - / ds — — — . (7) 
J 1 + u ■ e 2 

then xc = Xend+Xint- A similar result has been found 
by direct integration of eq. . We are rather remarkably 
back to a problem in statistical mechanics which is very 
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close that of the Fuller formulation, eq. (0, of the writhe 
fluctuations. We convert eq. Q) to spherical coordinates 
and find: 

f d(j) 

Xend= / d,S (1 - COsd) — (8) 

The singularity at 1 + u • e z = corresponds to the limit 
cos 8 = — 1 . 



E. Twisting 

Experiments are sensitive to the sum of the writhing 
and twisting fluctuations in a polymer. If we now allow 
the excitation of the torsional modes in the polymer of 
figure |21 the total rotation between the two ends is the 
sum of the writhing angle \c and rotation due to in- 
ternally excited twisting motions. In order to compare 
the twisting fluctuations to writhing fluctuations, we de- 
fine it = C/ksT the torsional length, where C is the 
torsional modulus 0. This twisting mode is unmodified 
by the writhing geometry, at least in the simplest models 
of DNA elasticity. The mean square twisting angle, \t 
is given by (x|) = j- t 



III. ARTEFACTS IN THE FULLER 
FORMULATION OF WRITHE 

Recent analytic work [|| Q has been based on the sim- 
plified formulation of eq. ®. It is much easier to treat 
analytically than the full double integral in equation (0) 
due to a direct mapping onto quantum mechanics: The 
bending energy of a stiff beam in the slender body ap- 
proximation of elasticity is given by 

A f fdt(s)\ 2 , 

where A is the bending modulus linked to the persistence 
length by £ p = A/k B T 

To calculate the partition function one must now sum 
over all paths 

Z = e- E ' kBT (10) 

paths 

The sum for the partition function is clearly closely re- 
lated to path integrals studied in quantum mechanics. 
Formally the energy in eq. 10 looks like the kinetic en- 
ergy of a free particle moving on a sphere. From the sum 
of paths in equation (|10|l one derives a Fokker-Planck 
equation which is entirely analogous to the Schrocdinger 
equation for a particle on a sphere: 



where V 2 is the Laplacian operator on the sphere. Here 
P(t, s) is the probability of finding the chain oriented in 
the direction t at the point s. As a function of s the 
vector t "diffuses" with diffusion coefficient which varies 
as l/£ p . 

Let us now consider the writhe calculated with Fuller's 
formula for such "diffusing" paths. There is a singularity 
in equation © near t = — e z . In the neighborhood of 
this direction the integral measures twice the winding of 
the random walk about the point — e z . The winding of 
a random walk has singular properties in the two dimen- 
sional plane and on a sphere jig. I n particular there is 
logarithmic divergence in the winding properties in the 
continuum limit. It is this winding number divergence 
that was picked up in the analytic calculation of [8j and 
necessitated an intermediate scale cut off in the calcula- 
tion. 

One of the principle conclusions of this paper is that 
these winding number singularities are not present in the 
distribution of the writhe determined experimentally: As 
shown above in the case of constrained linking number it 
is the Calugareanu- White expression for the writhe which 
gives the exact twisting angle. The Fuller expression dif- 
fers by an arbitrary factor of 4n7r due to the winding 
number singularities not present in the original formu- 
lation. We shall demonstrate numerically that it is the 
passage from © to ^ which introduces these singular- 
ities. Use of © will be shown to lead to a substantial 
error in the calculation of xc with an overestimate by a 
factor of 2.5 for a discretization corresponding to DNA. 



IV. NUMERICAL METHODS 

Given the difficulty of treating the full Calugareanu- 
White expression for the writhe analytically we decided 
to proceed by numerical exploration of the distributions 
of writhe implied by the Calugareanu- White and Fuller 
formalisms. In particular we look for the Cauchy tail 
predicted analytically |8( in the writhe distribution func- 
tion. The existence of such a tail implies the absence 
of a reasonable continuum limit of the wormlike chain 
and continuous evolution of the response functions as a 
function of a microscopic cut off. We shall conclude that 
the Calugareanu- White formulation remains finite even 
in the continuum limit. 

In our numerical investigations we shall be particularly 
interested in the fluctuations in the writhe of the open 
chain, as a function of the tension. Due to the algorithm 
used in generating the chains we are unable to generate 
chains in an ensemble with an imposed torsional couple. 
Our results are thus always for T = 0. The torsional 
fluctuation are, as usual, related to the linear response of 
the chain via the fluctuation-dissipation theorem. 
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A. Numerical Calculation of Writhe 

A number of methods are available for the calculation 
of the writhe of a discretized polymer 0] . We calculate 
the discretized versions of the integrals of eq. J2J an d 
equation J3) by recognizing that they are both areas on 
a unit sphere. In the Fuller formulation it is the area 
enclosed by the curve t(s). For the discretize chain the 
tangent curve becomes a series of link directions, t, , cor- 
responding to points on a sphere. One connects these 
points by geodesies and then sum over the area of the 
triangles formed by two successive tangent vectors and 
e z . 

XF =^A(e z ,t„,t n+ i) (12) 

n 

where A is the area of the spherical triangle defined by 
the three vectors. 

We have already noted that the full expression, eq. (0) 
corresponds to the area swept out by e(s,s'). For two 
links forming a discrete chain this defines a spherical rect- 
angle. We calculate its area by decomposing it into two 
spherical triangles. 

In both calculations we calculate the area of a spherical 
triangle by using PHuilier's expression 

_4 = 4arctan( ( tan((o + b + c)/4) x (13) 
tan((c - a + 6)/4) X 
tan((c + a — 6)/4) x 
tan((a + 6-c)/4)) 1/2 ) 

where, a, b and c are the lengths of the sides of a spherical 
triangle. For our purposes the triangle has to be oriented, 
the area can be either positive or negative. 

A useful cross check in the programming is that for 
any polymer the modulo relation of eq. (JIJ must be sat- 
isfied despite very different intermediate results in the 
calculation. Numerically we found that the equality held 
to within 10 -10 , when working in double precision when 
the extended Calugareanu- White definition of writhe was 
compared with the Fuller formulation. 

Since the double integral of equation is reduced to 
a double summation this step takes 0(N 2 ) operations 
for a chain of N links. For the long chains studied in our 
simulations it is by far the slowest step in the calculation. 



B. Link Sector Choice 

The generation of long unknotted chains is numeri- 
cally difficult; In our simulations we used an ensemble of 
chains with only the constraint Ck = which, as shown 
above, is the minimum constraint needed for unbounded 
torsional response. Our result can then be directly com- 
pared with existing analytic theories which do not impose 
any topological constraints. 



The ensemble of Ck — chains may contain knots. 
While this ensemble may appear physically "unreason- 
able" one must not forget that knots are rare in the chains 
that we shall study. It has been noted [2(| that for a flex- 
ible chain unknotted configurations dominate the statis- 
tics of chains even several hundred Kuhn lengths long. 
For the chain lengths that we work with in this paper 
the contamination coming from such knotted configura- 
tions should be weak. At the end of the paper we present 
a partial investigation of the influence of knots. We find 
that they do not modify our conclusions as to the nature 
of the continuum limit for writhing chains. The errors 
due to the use of a knotted ensemble are much smaller 
than the differences between the Fuller and Calugareanu- 
White formulation of the writhe. 



C. Chain generation 

In order to use the expressions for the writhe given 
above we are interested in chains in which the initial 
and final tangents are parallel (though the writhe for 
non parallel configurations does have a simple general- 
ization |l4f) but for which there is no constraint on the 
final position of the chain. Rather than using a conven- 
tional Monte-Carlo algorithm to generate chains we used 
a simple "growth algorithm" . 



1. Zero force 

We wish to grow chains of length L, persistence length 
£ p using a series of links of length b. In the absence of 
tension we generate chains by starting from a single link 
in the e z direction at the origin. We then successively 
add links to the chain with small random angle incre- 
ment ao ~ ybjtp to produce a single realization of an 
equilibrated semiflexible chain. It is almost certain that 
this chain does not satisfy the boundary conditions on 
the tangent thus we continue growing until the final tan- 
gent is parallel to the initial tangent to within an angle 
small compared with «o- We keep the polymer in our en- 
semble if the length of the polymer is less than 1.05 x L, 
otherwise the whole configuration is rejected and the pro- 
cess restarted from the first link. The configurations are 
then used to calculate the writhing distributions. The 
curves that we generate are somewhat "imperfect" since 
they are due to a mixture of lengths. This admixture of 
chain lengths plays no role, however, in our analysis of 
the asymptotic distribution of writhe. There is no self 
avoidance in this code; it can be shown from a Flory ar- 
gument that self avoidance is a weak effect in semiflexible 
chains of moderate length. 
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2. Finite force 

In the presence of an external force the algorithm is 
slightly more complicated. We proceed by noting that 
the partition function of a chain under tension can be 
expressed in a very similar manner to the partition func- 
tion of a flexible polymer in an external potential |2l| . 
We proceed by simulating the equation 

BZ 1 

_ = ^ + /cos(*)Z, (14) 

where 9 is the angle between the direction of the force 
and the local tangent to the polymer, and / = F/ksT. 
Z(s,9) corresponds the number of configurations in 
which the chain points in the direction 9 after a distance 
s. 

As proposed in [22] we introduce a pool of several chains 
which we grow simultaneously. As each link is added 
there is angular diffusion as described above and a sec- 
ond process of birth or death of chains in the pool in 
order to account for the force. If / cos(6>) is positive then 
it is considered to be a growth rate for reproduction of 
chains in the pool. If /cos(#) is negative the chain is 
stochastically destroyed with the appropriate probabil- 
ity. We also manage the total pool size as in |22j • At the 
end of a pool growth, we destroy the chains that do not 
satisfy the condition for the tangent vector to be parallel 
to within ao at both extremities. 

There are several sources of error possible with the al- 
gorithm. The most difficult to evaluate is the effect of the 
finite pool size. The result of a single run is an ensem- 
ble of several configurations together with a total weight 
coming from the management of the pool size. For suf- 
ficiently large pool sizes this weight is the same for each 
realization of the growth process. For small pool sizes, 
however, this weight undergoes important fluctuations. 
To calculate a correlation function from an ensemble of 
pools we chose to select a single chain from each pool and 
performed a simple average over at least 10, 000 pools. 
Since we had no, a priory method of estimating errors 
from this procedure we experimented with the pool size 
for several different values of the force. We found that 
even when varying the pools size from as low as 20 chains 
to 2000 chains the estimates of the mean square writhe 
were stable within a few per cent. In our production runs 
we chose a value of 50 chains per pool. An alternative 
procedure would weight each pool according to the true 
variation of the total weight of each simulation. 

A systematic difference between a discretized chain 
and a continuous curve also occurs. Our principle aim 
is to understand experiments on DNA0,0 therefore we 
choose £ p — 53 nm. We chose for b the half pitch of a 
single helix, 1.8 nm, yielding l v ~ 306. This is also com- 
parable to the diameter of the molecule. In order to study 
the convergence of the writhe to the continuum limit we 
shall also perform some simulations with £ p /b 30 (see 
section IVI Cll . 



V. ASYMPTOTICS OF WRITHE IN OPEN 
FLEXIBLE CHAINS 

Before presenting our numerical results on semiflcxi- 
ble chains we wish to explore the scaling behavior of the 
writhe of a polymer described as a freely jointed chain 
with N links. This allows us to understand the length 
scales and the structures which are important in deter- 
mining the writhe of a long molecule. Many of these 
results are already known for closed chains however we 
wish to demonstrate that the end corrections in eq. Q 
are subdominant in the limit of very long chains. 

A. Scaling arguments for flexible chains 

We shall start with the interpretation of the writhe as a 
signed area swept out by the vector e(s, s'): Consider two 
links i and j of length a separated by a distance Rij . This 
area scales as Ay ~ ±(a/Rij) 2 when R^ ^> a. When 
a <C Rij the area is bounded above by 2tt. Clearly when 
one averages over random walks the integral in equation 
gives zero. The mean squared writhe can be estimated 

as 

(xc> ~ £ A lnRa) (15) 

pairs 

where P(Rij) is the pair distribution function for the 
polymer. For a three dimensional Gaussian polymer this 
function scales as 

P(%) ~ 1/Rij (16) 

for lengths smaller than the radius of gyration, Rq ~ 
VaL of the polymer. Approximating the sum by an in- 
tegral we find the internal contribution to the writhe, 

This integral converges at large distances, but is diver- 
gent at small distances. We conclude that the writhe in- 
tegral is dominated by the cut off scale a. For a semiflcx- 
ible polymer the cut-off a corresponds to the persistence 
length l p . It is thus the structure at this length scale 
which dominates the writhe of the molecule. The aver- 
age crossing number is defined in a very similar manner 
to the writhe except the unsigned area |Ay| is summed 
over rather than the signed area. This has a non-zero 
mean and scales in the following way in a Gaussian poly- 
mer 

leading to a logarithmic divergence. Thus in the crossing 
properties of an arbitrary projection of the polymer we 
expect all length scales are important. 
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Finally, eq. J7J), the end correction to the writhe is a 
sum of random areas Aj ~ ±a/i?j where Ri is now the 
distance between the end of the polymer and the single 
link i. The corresponding estimate of the mean square 
writhe is thus 



(Xend) 



1 



a a 
~R2R 



R 2 dR ~ log ■ 



(19) 



giving a logarithmic contribution with the structure of 
the whole molecule being important. 

When self avoidance is introduced in the problem we 
use the result that P(R) ~ l/i? 4 / 3 ^ to show that the 
integral in equation l)17|l is still dominated by short length 
scales. The results for the end correction and average 
crossing number are however modified. They too be- 
come sensitive to structure at short wavelengths in the 
polymer. Thus from eq. i|17fl and eq. I|19|) the end correc- 
tions remain small for long chains and are subdominant, 

(Xint) > (Xlnd)- 

Given the importance of this end correction in the in- 
terpretation of the experiments we now present a more 
rigorous study of the end correction and confirm its sub- 
dominant nature compared with the internal contribu- 
tions to the writhe. 



B. Magnitude of End Corrections 

In order to further study the scaling behavior of the 
end corrections we examine the limit L/£ p large. We thus 
study the problem of a freely jointed chain rather than 
the semiflexible chain and disprove arguments of [24j that 
the dominant singularities in the writhe of a semiflexible 
chain come from ends due to the formal analogies be- 
tween eq. @ and eq. (JJJ. Numerical results (not shown 
here) on semiflexible chains lead to the same conclusions. 

Let us now make a hypothesis that the asymptotics of 
Xend are dominated by the largest contributions to the 
integrand of eq. (JHJ then 



Xend 



ds 



ds 



(20) 



which is the winding number of the polymer about e z . 
The winding number of a semiflexible random walk, W, 
about an infinite line scales as 



W 2 ~ log 2 (L/£ p ) 



(21) 



where L is the length of chain and £ p the persistence 
length. The Cauchy singularity appearing in this prob- 
lem is regularized by the stiffness of the chain, £ p . 

To check this hypothesis (and to improve on the very 
rough scaling analysis given above, eq. 1|19|0 we gener- 
ated, figure 01 a large number of random walks for a 
freely jointed chain. We then plot (x' 2 n d) as a frac- 
tion of logL and look for a straight line. We conclude 
that the end corrections to the writhe are comparable to 
(xtnd) = l°g 2 (-k/^p)- They are negligible compared with 
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FIG. 4: Variation of end correction, (x'i n d) i with L show- 
ing scaling identical to that of the winding of a random walk 
about an infinite line. Freely jointed chains with up to 320,000 
links. 



the internal contributions which vary as (xfnt) ~ i^/^p) 
for L > £ p . 

Despite the similarity between eq. J7J and eq. J3J we 
find very different scaling for (x%) ~ L\og(£ p /b) and 

(xlnd) ~ l°g 2 {L/£ p ). At first sight this might seem 
rather surprising, however in the Fuller formulation one 
averages over realizations of two dimensional random 
walks in the surface of the sphere whereas in equation 
@ we average over three dimensional random walks pro- 
jected onto a sphere. The statistical weights are different 
even if the functions are similar. 



VI. WRITHE DISTRIBUTION OF 
SEMIFLEXIBLE POLYMERS 

A. Short Molecules 

For short filaments of length L -C £ p the distribution of 
writhe calculated with the extended Calugareanu- White 
formula and the Fuller formula are indistinguishable; am- 
biguities due to winding about the pole are exponentially 
rare. The writhe distribution of a open polymer with par- 
allel tangents at each end in the limit L/£ p <C 1 is given 
by the Levy H| formula for the distribution of the 
area enclosed by random walk in a plane 



P(xc) 



f 



1 



2L cosh 2 (xcip/L) 



(22) 



By generating an ensemble of 10 6 short chains and bin- 
ning the writhe we verified that our code was able to 
reproduce this result. 
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B. Long molecules, zero tension 

We have characterized the evolution of the writhe prop- 
erties of a chain as a function of its length, L. We plot, 
figure |S] top, (xc) as a function of L for £ p fixed. For 
small L, eq. I|22|) . we have (xc) ~ an d f° r l° n g chains, 

we have (xh) ~ ^> ec l- CHI- 

As a second characteristic of the writhe dis- 
tribution, figure [3] bottom, we consider p4 = 
((x — (a;)) 4 ) / ((x — (s)) 2 ) related to the kurtosis, cal- 
culated for a probability distribution If p is Gaus- 
sian, then p4 = 3. For the distribution eq. I)22|l /?4 = 
21/5. We see numerically that p4 — 21/5 for small L 
and tends slowly to 3. However, there is a large peak 
around L/£ p — 2. The very strong non-monotonic be- 
havior in figure bottom is most striking; we shall now 
explain the origin of this feature. 

Let's call back-facing segments the sections of the chain 
along which t ■ e 2 < 0. We define, ng, the number of 
such segments for a chain. If we compute p^ with only 
configurations for which n$ — the peak disappears. We 
understand that the peak is due to the formation of loops 
in the chain which can, sometimes, completely dominate 
the writhing properties of a chain. Such important fea- 
ture in the distribution are clearly missed in any Monge 
description of a chain. 

From these simulations, we also estimated the coeffi- 
cient of proportionality K between (xc) an( i L/£ p for 
very long chains when the distribution of writhe has con- 
verged to a Gaussian: 



En 



K = 2.85 ± 0.03 



(23) 



the coefficient is similar in magnitude to that found 
|27| for closed chains, even though our model for the chain 
is somewhat different. 
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FIG. 5: Top: (xc) as a function of L. Bottom: pi as a 
function of L/£ p . When L/£ p — » we find pi ~ 21/5. There 
is a peak around L/£ p ~ 2, due to the presence of a small 
number of loops (typically one or two), which enlarge the 
probability distribution of the writhe. For larger values of 
L/lp p4 — * 3 corresponding to a Gaussian distribution. 
corresponds to chains with ns = 0. 100,000 chains for each 
point, £ p = 1000b for short chains and £ v = 506 for long 
chains. 



C. Convergence of the writhe distribution 

In order to characterize the asymptotic distribution of 
the writhe, and study the importance of the Cauchy tail 
in Fuller formulation, we have performed simulations on 
a series of chains of length L = 8£ p . With chains of this 
length knots remain rather rare whilst the energetic bar- 
rier needed for a chain oriented in the direction e 2 to wind 
about the direction — e z is only a few fc^T. We are thus 
sensitive to the winding singularities of the Fuller formu- 
lation. In our simulations we vary the discretization so 
that there are L/£ p = 10, 30, 100, 300, 900, 2700 links 
per persistence length and we generate 200000 indepen- 
dent configurations for each value of L/ell p . We plot the 
variance of xf and xc as a function of the discretization 
in figure [BJ We observe a continuous evolution showing a 
logarithmic divergence for the Fuller formulated writhe, 
and a convergence to a stable value for the Calugareanu- 
White form for quite moderate values of L/£ p . The ratio 



between (Xf) anc ^ (Xc) f° r DNA value of £ p = 306 is 
2.4. We also see that the difference between a chain with 
30 links per persistence length with the continuum limit 
is small (about 3%). 

A divergence of the angular fluctuations implies the 
breakdown of linear response in the continuum limit. 
Since, however, (xc) converges to a finite value we con- 
clude that a microscopic cutoff is not needed to render a 
torsionally stiff chain finite. 

D. Long tense molecules 

Until now we have ignored the effect of tension on the 
configuration of the DNA, except to remark that even 
very low tensions justify the use of a conserved linking 
number in the interpretation of the experiments. In this 
section we indicate how the writhe of DNA varies in the 
presence of external forces. 
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FIG. 6: Variance of the writhe computed with two differ- 
ent formulations as a function of discretization. \f evolves 
logarithmicly, whilst \c converges to a constant. DNA corre- 
sponds to e p /b ~ 30, for which ( X f) / (xc) = 2.4. L/l p = 8, 
200000 chains per point. The statistical error is smaller than 
symbols' size. 




FIG. 7: (Xf) Pp/L as a function of the scaled tension. At 
large forces the result converges to eq. Il2(it . At low tensions 
the writhe is independent of the force. Simulations for £ p = 
2506. Curves for L/£ p = 0.25, 0.5, 1, 2, 4, 8, 16, 32 (bottom 
to top). The straight line is eq. 1261 . 



When a molecule is under high tension the molecule is 
largely aligned parallel to the external force. We can use 
a simplified, quadratic form for the Hamiltonian [(| 



r(t_L At_L -e z ) 2 +k B Tft > ± \ ds 



E = \ J (9t± 

where tj_ describes the transverse fluctuations in the di- 
rection of the molecule. We can find the mean squared 
writhing angle using the usual methods of equilibrium 
statistical mechanics. 



(xl) = 



A short calculation gives 



d 2 log Z 
d {fiVf 




(25) 



(26) 



This expression can only be expected to be valid for forces 
such that ftp > 1. 

To estimate the writhe of a molecule at low forces we 
return to the remark above that the internal contribu- 
tion to the writhe is dominated by structure occurring 
at the scale £ p . Under low tensile forces the structure of 
a polymer is unchanged out to a length scale if = 1//- 
We conclude that under low tension, when £f 3> £ p the 
writhe of a molecule becomes independent of its degree 
of elongation and thus the tension, /. It is only un- 
der the highest forces when the semiflexible nature of 
the molecule is sampled that we see an evolution of the 
writhe with force. 

To validate our code for generating tense molecules we 
performed a series of simulations with a chain of per- 
sistence length £ p = 250&, results are shown on figure [7| 



using the Fuller expression for the writhe. At high forces 
we see that the curves converge towards the law eq. <|26[) 
and that at low forces the curves saturate as expected 
from the above scaling argument. Somewhat surprising 
however is the rapid crossover which occurs for forces 
comparable to / ~ l/£ p when L/£ p 3> 1. In the neigh- 
borhood of this force there is a pronounced "shoulder" on 
the curve. Convergence to the long chain limit is rather 
slow. It is not until L ~ 60£ p that we see a saturation 
in writhing curves. This slow convergence can be under- 
stood rather easily by noting that any chain within £ p 
of the surface of the polymer coil is in a region of lower 
than average density. 

We have performed a series of simulations on differ- 
ent levels of discretization of the polymers. We find that 
as the discretization becomes coarser the large shoulder 
dominates over the law in 1/s/f for the mean square 
writhe. This is illustrated in figure [S] where we use 30 
links per persistence length for the discrete chain. The 
chain of length L = £ p still displays a regime in agreement 
with eq. (|26[1 . However with longer chains the regime in 
l/yff is overwhelmed by the crossover to the low tension 
regime. Eq. (|26|1 substantially underestimates the writhe 
fluctuations in the domain out to f£ p = 10 correspond- 
ing to very large forces of ~ lpN. At such forces other 
corrections come into play, including the chiral nature of 
the DNA chain. 

It is interesting to note that the Fuller formulation 
gives results which are useful over a larger window of 
forces. The Fuller and Calugareanu- White curves are 
very similar down to tensions f£ p ~ 2; analytic calcu- 
lations based on the Fuller formulation can be expected 
to be useful for tensions larger than 0.2 pN. At very 
low forces the fluctuations are rather different. The ratio 
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FIG. 8: (xl) (p/L (left) and (xc) i P /L (right) as a function 
of the scaled tension. For longer chains the expected regime 
in \/\ff is hidden by a large shoulder from the crossover to 
the low force regime. Simulations for £ p = 306. Curves for 
L/l p = 1, 2, 4, 8, 16 and 32 (bottom to top). The straight 
bold line is equation 1261 . At zero force, the ratio (xf) I (Xc) 
tends to 2.5 as L/£ p grows. The horizontal dashed line is 
the amplitude of fluctuation due to twist fluctuations when 
(■t/ip = 1-5. The vertical dashed line corresponds to F = 
0.1 pN for DNA. 



(Xf) I (xc) i s about 2.5; use of the Fuller formulation 
strongly over estimates the torsional response functions. 



E. Origin of the shoulder 

To understand the origin of the shoulder which appears 
in figure 6 we have classified the different configurations 
as a function of the number of back-facing segments and 
computed (xc) f° r subsets with ns fixed. The results 
are shown in figure for £ p /b = 30. When f£ p ~ 1 we 
see that the curves for different values of ns separate. 
ns = curves do not exhibit a shoulder. 

The origin of the shoulder in figure |H1 is essentially the 
same as the origin of the peak in figure Under a large 
force, the configurations with back-facing segments are 
statistically rare because they have a large energy cost. 
Their contribution to (xc) i s > however, important. 



VII. WRITHE AND KNOTS 

Clearly, self avoidance of the chain implies that topo- 
logical invariants are constrained |24|. including Ck. An- 
other topological invariant that should be constrained is 
the knot configuration. The algorithm that we have used 



O O o 



o O 









I P /L 




I P /L 


>< <x c 2 V 




I p/L 


+ <*c\=2 




FIG. 9: (xc)£p/L as a function of f£ p . Circles all chains. 
Triangles, averaged from no back-facing segment. We have 
plotted the curves for ns = 1, X and ns = 2, (+). Each point 
computed from 300,000 configurations of length 8£ p with £ p = 
306. Equation solid line 



until now leads to an ensemble containing both knotted 
and unknotted chains, even if Ck — 0. We now present 
a preliminary investigation on the influence of knots on 
the distribution of writhe. This question is related to the 
problem of the closure of the chain: If the polymer is 
allowed to pass through the line to infinity, a knot may 
appear, change or disappear. We thus continue with the 
assumption that the tension of the DNA is sufficiently 
high that the extended construct of figure 2 remains in 
the same state throughout an experiment. 

In order to select the knotted chains, we have com- 
puted the Jones polynomial V [2^, [2£|. The choice of 
the Jones polynomial is justified by the Jones conjecture, 
that is Vk = 1 if and only if K is the trivial knot. We 
have used the algorithm of Kauffman p^.l30| to calculate 
Vk and perform the classification of chains in function of 
their knots. Recently it has been noted that probabil- 
ity of knot formation is small [2(J in short chains. Our 
simulations confirm this point : for chains of 8 persis- 
tence lengths the proportion of knotted chain is around 
(5 ± 1) x 10~ 4 when / = 0. 

To interpret a micromanipulation experiment one 
should perform averages over an ensemble of chains with 
the same given knot. In practice one hopes the experi- 
ment is performed with the trivial knot. The probabil- 
ity distribution of the writhe angle xc has to be mod- 
ified, because we do not count the knotted configura- 
tions. We have computed (xc)tk> wnere the label "tk" 
stands for "trivial knot" and compared it to (xc)j f°r 
different lengths with no force. Results are plotted on 
figure E| The knot configurations that our algorithm 
generated lead to an overestimated value for (xc) ■ The 
correction is about —7.5% for the longest chains that 
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FIG. 10: Fluctuations of \c f° r all chains, compared with 
unknotted chains, as a function of L/£ p . T (xc) tk i s be- 
low (xc) ■ The knot correction lowers the fluctuations. For 
L/£p = 80 the correction computed is —7.5%. (3<r) shown for 
(Xc) tk ; error for (xc) sma ll than symbol size. 



strand so that under torque we expect an series of equiv- 
alent low energy states separated by a potential barrier. 
For the usual bead sizes and forces used in the experi- 
ments this barrier is very high and the extended linking 
number is conserved. This allows the use of an extended 
Calugareanu- White formalism in the calculation of the 
bead rotations. We find that the Fuller and CaLugareanu- 
White formulae give substantially different distribution 
functions for the torsional fluctuations due to writhe. In 
contrast to Bouchiat et al. we find that the topologi- 
cally confined DNA chain does not need an intermediate 
scale cut off to render the response functions finite. The 
mathematical problems as to existence of the torsional 
response functions occur at long wavelengths; imposing 
a short wavelength cut off is the wrong solution to this 
problem. 



we studied of L = 80£ p . The correction increases the 
disagreement between calculations based on Fuller's for- 
mulation of the writhe and the experimental curves. It 
is particularly instructive to compare figure 1101 and fig- 
ure for the case L/£ p = 8. Removing knots from the 
ensemble of chains has a negligible effect on the writhe 
for such short chains, however even for the most coarsely 
discretized chains with £ p /b = 10 the ratio (Xf)/(Xc) 
is larger than 2. We have also studied the evolution of 
this correction with the force. The correction decreases 
rapidly with the force and is found to be zero for forces 
Ftp larger than a few ksT. 

VIII. CONCLUSIONS 

In this paper we have shown that the standard ex- 
perimental geometry does not torsionally confine a DNA 



In experiments we expect several distinct regimes when 
working with beads of size d t p . For very low forces, 
F < ksT/d torsional fluctuations are unbounded and it 
is not possible to define the torsion-force-extension char- 
acteristics. In the regime ksT/d < F < ksT/tp torsional 
fluctuations are bounded but must be calculated using 
the full double integral representation of the writhe. The 
Fuller formulation, even with an additional cut off, sub- 
stantially overestimates the writhe contribution to the 
torsional response. When F£ p ~ ksT a simple theory 
based on a Monge representation expanded to quadratic 
order is unable to fit the data due to strong corrections 
to scaling; here an expression based on the Fuller for- 
mula can be expected to give a better description of the 
response. Finally, for very large forces, F ^> ksT/tp, the 
Calugareanu- White and Fuller formulae give the same 
result, however other effects which are neglected here be- 
come important; a full theory must treat the chiral nature 
of DNA and force induced denaturation. 
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